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Calculational results are presented for the fine-structure splitting of the 2^P state of helium and 
helium-like ions with the nuclear charge Z up to 10. Theoretical predictions are in agreement with 
the latest experimental results for the helium fine-structure intervals as well as with the most of 
the experimental data available for light helium-like ions. Comparing the theoretical value of the 
interval in helium with the experimental result [T. Zelevinsky et al. Phys. Rev. Lett. 
95, 203001 (2005)], we determine the value of the fine-structure constant a with an accuracy of 31 
parts per billion. 
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I. INTRODUCTION 



The fine structure splitting of the 2^P state in helium plays a special role in atomic spectroscopy because it can be 
used for an accurate determination of the fine structure constant a. This fact was first pointed out by Schwartz in 
1964 [1]. The attractive features of the fine structure splitting in helium as compared to other atomic transitions are, 
first, the long lifetime of the metastable 2^Pj levels (roughly two orders of magnitude larger than that of the 2p state 
in hydrogen) and, second, the relative simplicity of the theory of the fine structure. Schwartz's suggestion stimulated 
a sequence of calculations [2-5] , which resulted in a theoretical description of the helium fine structure complete up 
to order ma^ (or Ry) and a value of a accurate to 0.9 ppm [6]. 

The present experimental precision for the fine-structure intervals in helium is sufficient for a determination of a 
with an accuracy of 14 ppb [7, 8]. In order to match this level of accuracy in the theoretical description of the fine 
structure, the complete calculation of the next-order, ma^ contribution and an estimation of the higher-order effects 
is needed. The work towards this end started in 1990s and extended over two decades [9-19]. In 2006, the first 
complete evaluation of the ma^ correction to the helium fine structure was reported [20]. However, the numerical 
results presented there disagreed with the experimental values by more than 10 standard deviations. 

In our recent investigations [21, 22], we recalculated all effects up to order ma^ to the fine structure of helium 
and performed calculations for helium-like ions with nuclear charges Z up to 10. The calculations were extensively 
checked by studying the hydrogenic {Z oo) limit of individual corrections and by comparing them with the results 
known from the hydrogen theory. We found several problems in previous studies. As a result, the present theoretical 
predictions are in agreement with the latest experimental data for the fine-structure intervals in helium, as well as with 
the most of experimental data available for light helium-like ions. Comparison of our theoretical prediction for the 
2'^Po — 2'^Pi interval in helium (accurate to 57 ppb) with the experimental value [7] (accurate to 24 ppb) determines 
the value of the fine structure constant a with an accuracy of 31 ppb. 

The calculation of the ma^ correction for the fine-structure splitting of light helium-like atoms was reported in our 
recent Letter [22]. In this paper, we present an extended description of the ma^ correction and a detailed term-by- 
term comparison of our results with independent calculations by Drake [18] for helium and by Zhang et al. [12] for 
helium-like ions. 



II. THE SPIN-DEPENDENT mcx' CORRECTION 



The ma^ correction to the fine-structure splitting of a two-electron atom can be conveniently separated into four 
parts. 
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The first term in the brackets above combines all terms with InZ and In a [10-12, 15, 20], 
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where f = fi — f2, Hq and Eq are the Schrodinger Hamiltonian and its eigenvalue, and H^^^ is the spin-dependent 
part of the Breit-Pauli Hamiltonian, 
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The second part of i''-^-' is induced by effective Hamiltonians to order ma^ . They were derived by one of us (K.P.) in 
Refs. [20, 21]. (The previous derivation of this correction by Zhang [10, 11] turned out to be not entirely consistent.) 
The result is 
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The Hamiltonian Hq is induced by the two-photon exchange between the electrons, the electron self-energy and the 
vacuum polarization. It is given by [20] 
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Here, the terms with In Z compensate the logarithmic dependence implicitly present in expectation values of singular 
operators 1/r^ and 1/r^, so that matrix elements of Hq do not have any logarithms in their 1/Z expansion. The 
singular operators are defined though their integrals with the arbitrary smooth function /, 
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where 7 is the Euler constant. 
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The effective Hamiltonian Hh represents tlie anomalous magnetic moment (amm) correction to the Douglas-Kroll 

ma^ operators and is given by [20] 
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The Hamiltonian -fffj'' is the ma^ amm correction to the Breit-Pauli Hamiltonian, 
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(ciC2 + C3) , 
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where ci = 1/2, C2 — —0.328 478 965 and C3 = 1.181241456 are the expansion coefficients of the free-electron amm in 
powers of (a/vr). 

The third part of is given by the second-order matrix elements of the form [20] 
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where H^^^ is the effective Hamiltonian responsible for the nonlogarithmic ma^ correction to the energy. 
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-^nfs spin-independent part of the Breit-Pauli Hamiltonian (with the term oc 5^{r) omitted since it does not 

contribute in our case), 
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and Hf^^ is the ma^ amm correction to H^^^ , 
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The fourth part of f*-^' is the contribution induced by the emission and reabsorption of virtual photons of low 

(7) 

energy. It is denoted as ' and interpreted as the relativistic correction to the Bethe logarithm. The expression for 
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TABLE I: Summary of individual contributions to the fine-structure intervals in helium, in kHz. The parameters [23] are 
a"^ = 137.035 999 679(94), cR-y^ = 3 289 841960 361(22) kHz, and m/M = 1.370 933 555 70 x 10"*. The values by Drake are 
taken from Table 3 of Ref. [18]. The label (+m/M) indicates that the corresponding entry comprises both the non-recoil and 
recoil contributions of the specified order in a. The uncertainty due to the value of a is not shown. 



Term i>oi vii Ref. 



ma*{+m/M) 29 563 765.45 


2 320 241.43 






29 563 765.23'' 


2 320 241. 42"^ 


[18] 


ma^(+m/M) 


54 704.04 


-22 545.00 






54 704.04 


-22 545.01 


[18] 


mo? 


-1607.52(2) 


-6 506.43 






-1 607.61(4) 


-6 506.45(7) 


[18] 


ma^m/M 


-9.96 


Q 1 




-10.37(5) 


9.80(11) 


[18] 


ma' \og{Za) 


81.43 


-5.87 






81.42*" 


-5.87*' 


[18] 


ma' , nlog 


18.86 


-14.38 




ma^ 


±1.7 


±1.7 




Total theory 


29 616 952.29 ± 1.7 


2 291 178.91 ±1.7 




Experiment 


29 616 951.66(70)= 

29 616 952.7(10)'* 
29 616 950.9(9)"= 


2 291 177.53(35)-^ 

2 291 175.59(51)^= 
2 291 175.9(10)3 





the original result was scaled to the present value of a. 
the original result was altered by the substitution ln(a) — >• \n{Za) in the terms proportional to ln(a), in order to comply 

with the present result for the logarithmic mci' contribution. 
" Ref. [7]. Ref. [24]. " Ref. [25]. ^ Ref. [8]. " Ref. [26]. 
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where 6 {■■ ■) denotes the first-order perturbation of the matrix element (. . .) by Hf^\ implying perturbations of the 
reference-state wave function, the reference-state energy, and the electron Hamiltonian. 



III. RESULTS: HELIUM 



The summary of individual contributions to the fine-structure intervals of helium is given in Table I. Numerical 
results are presented for the large uoi and the small V12 intervals, defined by 

1^01 = [Ei^'Po) - E{2''Pi)]/h, (15) 
= [E{2''P^) - EifP^)] /h . (16) 

We note that the style of breaking the total result into separate entries used in Table I differs from that used in the 
summary tables of the previous papers by K.P. et al. [20, 21]. Particulary, the lower-order terms listed in Table III of 
Ref. [20] and in Table II of Ref. [21] contained contributions of higher orders, whereas in the present work, the entries 
in Table I contain only the contributions of the order specified. 

A term-by-term comparison with the independent calculation by Drake [18] is made whenever possible. We observe 
good agreement between the two calculations for the lower-order terms, namely, for the ma'*, m.a^ , and ma^ correc- 
tions. However, for the recoil correction to order ma®, our results differ from Drake's ones by about 0.5 kHz for both 
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TABLE II: Individual contributions to the fine-structure intervals of helium- like atoms, in MHz/Z*. The label {+m/M) 
indicates that the corresponding entry comprises both the non-rccoil and recoil contributions of the specified order in a. For 
Z — 3, 7, and 10, a term-by-term comparison is made with the previous calculation by Zhang et. al. [12]. The results of 
Ref. [12] for the leading nia^ correction were scaled to the present value of a. The deviation for the ma'' (log) contribution is 
due to the difference in the expressions for this term. 



Z ma^{+m/M) ma^{+m/M) moL^ mo^mjM ?nQ^(log) ma^ (nlog) Total Ref. 



2 


1847.735 34 


3.419 00 


-0.10047 


-0.000 62 


0.005 09 


0.001 18 


1851.059 52(11) 


3 


1917.793 96 


3.249 78 


1.230 26 


-0.002 43 


-0.010 76 


0.018 01 


1922.278 81 (59) 




1917.793 97 


3.249 78 


1.230 25 




-0.010 27 




1922.26 (2) 


4 


1346.965 34 


1.943 84 


4.566 98 


-0.006 70 


-0.028 43 


0.046 48 


1353.4875 (39) 


5 


765.885 57 


0.685 51 


10.37447 


-0.01417 


-0.04139 


0.086 28 


776.976 (14) 


6 


270.38772 


-0.36757 


19.26647 


-0.02789 


-0.048 86 


0.139 52 


289.349 (37) 


7 


-139.085 57 


-1.229 55 


31.908 79 


-0.045 30 


-0.051 10 


0.209 03 


-108.294 (83) 




-139.085 58 


-1.229 55 


31.908 82 




-0.046 33 




-108.5(3) 


8 


-477.534 46 


-1.93791 


48.988 80 


-0.068 79 


-0.048 55 


0.297 85 


-430.30 (17) 


9 


-759.770 39 


-2.526 32 


71.203 90 


-0.093 96 


-0.04163 


0.409 16 


-690.82 (31) 


10 


-997.723 26 


-3.021 03 


99.25705 


-0.13744 


-0.030 76 


0.546 19 


-901.11(53) 




-997.723 25 


-3.021 03 


99.25705 




-0.021 29 




-901.5 



2 


1992.750 43 


2.009 94 


-0.50712 


-0.000 05 


0.004 72 


0.000 28 


1994.258 20(11) 


3 


1150.274 90 


-0.942 85 


-0.864 60 


-0.000 05 


-0.02216 


0.014 83 


1148.460 07(41) 




1150.274 91 


-0.942 85 


-0.86460 




-0.023 48 




1148.44 (2) 


4 


-384.659 15 


-4.448 24 


-1.389 63 


-0.000 06 


-0.045 39 


0.032 04 


-390.5104 (12) 


5 


-1739.328 53 


-7.320 66 


-2.393 83 


-0.000 04 


-0.05446 


0.046 61 


-1749.0509 (32) 


6 


-2838.550 28 


-9.580 33 


-3.994 54 


0.000 01 


-0.048 68 


0.056 88 


-2852.1169(77) 


7 


-3724.421 92 


-11.370 60 


-6.245 32 


0.000 08 


-0.029 03 


0.062 15 


-3742.005 (16) 




-3724.421 93 


-11.370 60 


-6.263 64 




-0.04190 




-3742.1 (3) 


8 


-4445.632 74 


-12.812 45 


-9.17416 


0.000 17 


0.003 27 


0.062 07 


-4467.554 (31) 


9 


-5041.009 23 


-13.993 89 


-12.79723 


0.000 25 


0.04705 


0.056 47 


-5067.697(55) 


10 


-5539.338 27 


-14.977 37 


-17.124 41 


0.000 38 


0.10127 


0.045 23 


-5571.293 (91) 




-5539.338 27 


-14.977 38 


-17.145 16 




0.075 67 




-5571.4 



[12] 



[12] 



[12] 



[12] 



[12] 



[12] 



intervals. The reason for this disagreement seems to be different for the large and the small intervals. For the large 
interval, the deviation is due to the recoil operator part, whereas for the small interval, it is mainly due to the mass 
polarization part (see discussion in Ref. [21]). 

Our estimates of the uncalculated higher-ordcir cfFcicts for helium are much larger than those in the previous studies 
[17, 18]. The previous estimates amounted to significantly less than 1 kHz for both intervals and were based on 
some logarithmic contributions to order mo? that were identified by analogy with the hydrogen fine structure. We 
now believe that the dominant mo? contribution might be of rclativistic origin. Our estimates of ±1.7 kHz for both 
intervals were obtained by multiplying the mo? contribution for the sum of vq\ -|- vyi by the factor of (Zoi)^ . 

Our result for the vq\ interval of helium agrees well with the experimental values [7, 24, 25] . For the vy^ interval, 
our theory is by about la away from the values obtained in Refs. [7, 26] but in agreement with the latest measurement 
by Hessels and coworkers [8]. Assuming the validity of the theory, we combine the theoretical prediction for the fgi 
interval in helium with the experimental result [7] and obtain the following value of the fine structure constant, 

(He) = 137.036 001 1 (39)theo (16)exp , (17) 
which is accurate to 31 ppb and agrees with the more precise results of Refs. [27-29]. 



IV. RESULTS: HELIUM-LIKE IONS 



Table H gives the summary of individual contributions to the fine-structure intervals of helium-like atoms with the 
nuclear charge number Z up to 10. We choose to present results for the intervals i/qi and z/o2 = z^oi + and not for 

vqx and v-y2, as is customary. The reason to consider 1^02 is that this interval is free from effects of the 2^Pi — 2^P-y 
mixing, which strongly affect the vq\ and vyi intervals. As a result of the absence of the mixing effects, all corrections 



6 



to 1^02 starting with the order of ma^ demonstrate a weaker Z dependence as compared to those to vqi and 1^12- 
The most drastic difference occurs for the a^w? jM correction: for Z = 10, this correction for ^02 is by 3 orders of 
magnitude smaller than that for vq\. 

The uncertainty of the theoretical values specified in Table II is solely due to uncalculatcd higher-order effects. Its 
estimation for helium was already discussed. For helium-like ions, we obtain the uncertainty by multiplying the raa^ 
contribution for the corresponding interval by the factor of (Zaf' . So, our error estimates are typically by a factor of 
XjZ smaller for the vq2 interval than for the 1^01 (or, equivalently, v\2) interval. 

For Z = 3, 7, and 10, Table II presents a term-by-term comparison with the previous calculation by Zhang et. 
al. [12]. We observe excellent agreement for the ma** and ma^ corrections. For the ma^ correction, the agreement is 
excellent in all cases except for the ^02 interval and Z = 7 and 10, where a small deviation is present. The results of 
the two calculations for the TOa^(log) correction are different, but this is explained by the difference in the expressions 
for this term. If we use the same expression as in Ref. [12], excellent agreement is found again. 

In Fig. 1 we plot our numerical results for the ma^ correction as a function of the nuclear charge number Z, together 
with the fit of the 1/Z expansion and with the asymptotical high-Z limit of this correction. The form of the l/Z 
expansion and the values of the first coefficient(s) are known. For the 1^02 interval, the leading term scales as Z^ and 
is calculated for hydrogen in Ref. [30]. For the 1^12 interval, there are additional and Z^ contributions due to the 
triplet-singlet mixing, which are obtained by expanding the following expression in 1/Z, 

|(2iPi|4^)|23Pi)|' 
'^^'"'^-i?o(23Fi)-i?o(2iPi)- 
The resulting asymptotic form of the nonlogarithmic ma^ correction is 

£:(^'"'°s)(j,g^)/^7 = 0.004045 - 0.015524/Z+O(l/Z2), (18) 
£:(^'"'°s)(j,g2)/zs = -0.021706 + 0(1/Z). (19) 

By fitting the 1/Z expansion of our numerical data for f ("^-"'og), we were able to reproduce well the values of the 
coefficients given above, which served as an important check of our calculations. 



0.8 




Z Z 



FIG. 1: Nonlogarithmic ma^ correction to the fine-structure intervals of helium-like atoms, for the uqi interval (left) and for the 
f02 interval (right). Black dots denote the numerical results, solid line stands for the fit of the 1/Z expansion, and the dashed 
(red) line indicates the asymptotic high-Z results. Note that the results in difi'erent graphs are scaled by difi'erent factors. It 
is Z'^ for the voi interval and Z^, for 1/02- The difi'erent Z scaling is the consequence of the triplet-singlet mixing effects. 

The comparison of the present theoretical predictions with experiment data is summarized in Table III. The 
agreement is very good in most cases. The only significant discrepancy is for Be^+, where the difference of 1.7 
standard deviations (a) is observed for 1/12 and that of 3.5 cr, for 1^02- It is important that for both the vqi and 
1^12 intervals there are experimental results available for helium-like ions, whose accuracy significantly exceeds the 
theoretical errors. These are the measurement of t'oi in helium-like nitrogen by Thompson at al. [31] and that of 
1^12 in helium-like fluorine by Myers et al. [32]. Good agreement of the present theory with these experimental 
results suggests that the theoretical errors (i.e., the uncalculated higher-order effects) were reasonably estimated. It is 
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TABLE III: Comparison of theoretical and experimental results for the fine-structure intervals of helium-like ions. Units are 
MHz for Li"*" and cm"'^ for other atoms. 



z 


Present theory 


Experiment 


Ref. 


vox 








3 


155 704.584(48) 


155 704.27(66) 


[33] 


4 


11.557 756(33) 


11.5586(5) 


[34] 


5 


16.198 21(29) 


16.203(18) 


[35] 


7 


-8.6731(67) 


-8.670 7(7) 


[31] 


8 


-58.791 (23) 


-59.2 (1.1) 


[36] 


10 


-300.58(18) 


-300.7(2.2) 


[36] 










9 


-957.886(79) 


-957.873 0(12) 


[32] 










3 


93025.266(34) 


93 025.86(61) 


[33] 


4 


-3.334663(10) 


-3.336 4(5) 


[34] 


5 


-36.463 787(66) 


-36.457(16) 


[35] 


8 


-610.392 3(42) 


-611.3(7) 


[36] 


10 


-1858.383 (30) 


-1858.3(1.5) 


[36] 



unfortunate that there are no experimental results with comparable accuracy available for the vqi interval. Since vqi 

is not affected by the triplet-singlet mixing effects, accurate experimental results for this interval in light helium-like 
ions would yield an improved estimate for the uncalculated higher-order effects in helium, thus increasing accuracy of 
determination of a from the helium fine structure. 

Comparing theoretical and experimental results for the fine structure of helium-like ions, one should keep in mind 
that the present calculation is carried out for a spinless nucleus, whereas the experimental results listed in Table III 
were performed for non-zero nuclear spin isotopes. For a nucleus with spin, the hyperfine splitting can usually be 
evaluated separately and employed for an experimental determination of the fine structure. This procedure, however, 
ignores the mixing between the hyperfine and the fine splittings. So, more accurate calculations should account for 
both effects simultaneously. 

In summary, the theory of the fine structure of helium and light helium-like ions is now complete up to orders moJ 
and o^rv? jM. Theoretical predictions agree with the latest experimental results for helium, as well as with most 
of the experimental data for light helium-like ions. A combination of the theoretical and experimental results for 
the 2^Po — 2^Pi interval in helium yields an independent determination of the fine structure constant a accurate to 
31 ppb. 

Support by NIST through Precision Measurement Grant PMG 60NANB7D6153 and by the Helmholtz Gemeinschaft 
(Nachwuchsgruppe VII-NG-421) is gratefully acknowledged. 



C. Schwartz, Pfiys. Rev. 134, A1181 (1964). 

M. Douglas and N. KroU, Ann. Phys. (NY) 82, 89 (1974). 

L. Hambro, Phys. Rev. A 5, 2027 (1972). 

L. Hambro, Phys. Rev. A 6, 865 (1972). 

L. Hambro, Phys. Rev. A 7, 479 (1973). 

M. L. Lewis and P. H. Scrafino, Phys. Rev. A 18, 867 (1978). 

T. Zclevinsky, D. Farkas, and G. Gabriclsc, Phys. Rev. Lett. 95, 203001 (2005). 

J. S. Borbcly, M. C. George. L. D. Lomhardi, M. Weel, D. W. Fitzakerley, and E. A. Hessels, Phys. Rev. A 79, 0605030(R) 
(2009). 

Z.-C. Yan and G. W. F. Drake, Phys. Rev. Lett. 74, 4791 (1995). 
T. Zhang, Phys. Rev. A 54, 1252 (1996). 
T. Zhang, Phys. Rev. A 53, 3896 (1996). 

T. Zhang, Z.-C. Yan, and G. W. F. Drake, Phys. Rev. Lett. 77, 1715 (1996). 
T. Zhang and G. W. F. Drake, Phys. Rev. A 54, 4882 (1996). 
T. Zhang, Phys. Rev. A 56, 270 (1997). 
K. Pachucki, J. Phys. B 32, 137 (1999). 



8 



K. Pachucki and J. Sapirstein, J. Phys. B 33, 5297 (2000). 
K. Pachucki and J. Sapirstein, J. Phys. B 35, 1783 (2002). 
G. W. F. Drake, Can. J. Phys. 80, 1195 (2002). 
K. Pachucki and J. Sapirstein, J. Phys. B 36, 803 (2003). 
K. Pachucki, Phys. Rev. Lett. 97, 013002 (2006). 

K. Pachucki and V. A. Yerokhin, Phys. Rev. A 79, 062516 (2009) [ibid. 80, 019902(E) (2009); ibid. 81, 039903(E) (2010)]. 

K. Pachucki and V. A. Yerokhin, Phys. Rev. Lett. 104, 070403 (2010). 

P. J. Mohr, B. N. Taylor, and D. B. NcwolL Rev. Mod. Phys. 80, 633 (2008). 

G. Giusfredi, P. C. Pastor, P. D. Natale, D. Mazzotti, C. do Mauro, L. Fallani, G. Hagel, V. Krachmalnicoff, and M. Inguscio, 
Can. .J. Phys. 83, 301 (2005). 

M. C. George, L. D. Lombardi, and E. A. Hessels, Phys. Rev. Lett. 87, 173002 (2001). 
J. Castillega, D. Livingston, A. Sanders, and D. Shiner, Phys. Rev. Lett. 84, 4321 (2000). 

D. Hanneke, S. Fogwell, and G. Gabrielse, Phys. Rev. Lett. 100, 120801 (2008). 

M. Cadoret, E. de Mirandes, P. Clade, S. Guellati-Khehfa, C. Schwob, F. Nez, L. Juhen, and F. Biraben, Phys. Rev. Lett. 
101, 230801 (2008). 

A.-M. Jeffrey, R. E. Elmquist, L. H. Lee, and R. F. Dziuba, IEEE Trans. Inst. Meas. 46, 264 (1997). 

U. Jentschura and K. Pachucki, Phys. Rev. A 54, 1853 (1996). 

J. K. Thompson, D. J. H. Howie, and E. G. Myers, Phys. Rev. A 57, 180 (1998). 

E. G. Myers, H. S. Margolis, J. K. Thompson, M. A. Farmer, J. D. Silver, and M. R. Tarbutt, Phys. Rev. Lett. 82, 4200 
(1999). 

E. Riis, A. G. Sinclair, O. Poulsen, G. W. F. Drake, W. R. C. Rowley, and A. P. Levick, Phys. Rev. A 49, 207 (1994). 
T. J. SchoU, R. Cameron, S. D. Rosner, L. Zhang, R. A. Holt, C. J. Sansonetti, and J. D. Gillaspy, Phys. Rev. Lett. 71, 
2188 (1993). 

T. P. Dinneen, N. Berrah-Mansour, H. G. Berry, L. Young, and R. C. Pardo, Phys. Rev. Lett. 66, 2859 (1991). 
N. J. Peacock, M. F. Stamp, and J. D. Silver, Phys. Scr. T8, 10 (1984). 



